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Abstract 

We examine the possibility of electroweak baryogenesis and dark matter in the nMSSM, a 
minimal extension of the MSSM with a singlet field. This extension avoids the usual domain 
wall problem of the NMSSM, and also appears as the low energy theory in models of dynamical 
electroweak symmetry breaking with a so-called fat-Higgs boson. We demonstrate that a strong, 
first order electroweak phase transition, necessary for electroweak baryogenesis, may arise in 
regions of parameter space where the lightest neutralino provides an acceptable dark matter 
candidate. We investigate the parameter space in which these two properties are fulfilled and 
discuss the resulting phenomenology. In particular, we show that there are always two light CP- 
even and one light CP-odd Higgs bosons with masses smaller than about 250 GeV. Moreover, 
in order to obtain a realistic relic density, the lightest neutralino mass tends to be smaller than 
Mz/2, in which case the lightest Higgs boson decays predominantly into neutralinos. 



1 Introduction 



The Standard Model (SM) provides an excellent description of all elementary particle interactions 
up to energies of order 100 GeV. However, there are several reasons to expect that an extension of 
the Standard Model description is needed at energies slightly above this scale. The most important 
of these is the precise mechanism for the generation of elementary particle masses. In the Standard 
Model, the Higgs mechanism provides a self-consistent parametrization of this mechanism, but relies 
on mass parameters that are sensitive to ultraviolet scales at the quantum level. Two ways to avoid 
this problem are to invoke super symmetry, or to break the electroweak symmetry dynamically. In 
both cases, the new physics provides an effective energy cutoff for the SM description that is of 
order of the weak scale. While the construction of supersymmetric extensions consistent with low 
energy observables is a relatively simple task 0121) this is much more difficult to do in models that 
break the electroweak symmetry dynamically due to the non-perturbative nature of the interactions 
in the symmetry breaking sector [3]. 

The Standard Model description also fails to answer two of the most important questions at 
the interface of particle physics and cosmology: the nature of the dark matter and the origin 
of the baryon asymmetry. A consistent dark matter density may easily arise in the presence of 
stable, neutral, weakly interacting particles with masses of order of the weak scale. In this sense, 
the observation of a non-vanishing dark matter density provides an additional motivation for the 
existence of new physics at the weak scale. The lightest supersymmetric particle (LSP) in models of 
low-energy supersymmetry with R-Parity conservation provides a natural source for the observed 
dark-matter density ^ |2] . 

There are several ways of explaining the origin of the matter-antimatter asymmetry. Some 
of them rely on physics at scales much larger than the weak scale. Leptogenesis, for instance, 
provides a mechanism for the generation of a lepton asymmetry via the decay of heavy Major ana 
neutrinos which explains the smallness of the neutrino masses in the See-saw framework. This lepton 
asymmetry is subsequently transformed into a baryon asymmetry via anomalous weak interaction 
processes: baryogenesis proceeds from leptogenesis [UIE]- Due to the high-energy nature of the 
fundamental interactions leading to leptogenesis, it is difficult to test the realization of this scenario 
by weak-scale experiments. 

Electroweak baryogenesis P2 provides an alternative to models of baryogenesis at high-energies, 
relying only on physics at the weak scale. The realization of this scenario demands that the 
anomalous baryon-number violating processes are out of equilibrium in the broken phase at the 
critical temperature of the electroweak symmetry breaking phase transition. This is possible only 
if the phase transition is strongly first order, or equivalently, 



where ip(T c ) is the Higgs vacuum expectation value in the broken phase at the critical temperature 
T c . Such a strongly first order transition cannot be achieved within the SM framework. 

Supersymmetric theories introduce new physics at the weak scale which may lead to a sufficiently 
strong first order phase transition. This can be achieved even within the minimal supersymmetric 
extension of the SM (MSSM). However, in the minimal case, this demands Higgs masses only 
slightly above the present experimental bounds and a light stop with a mass below that of the top- 
quark 7 . Therefore, this scenario is highly constrained and may be soon tested by the Tevatron 
collider experiments jSj- 
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As we will explain in the next section, there are good reasons to go beyond the minimal super- 
symmetric framework. In this work we shall study a minimal extension of the MSSM which includes 
a singlet field having restricted interactions. In Section |21 after motivating this extension, we will 
describe its most relevant properties. In Section |3] we will define the model precisely, set out our 
notation, and study the particle spectrum at zero-temperature as well as the collider constraints. 
Section 0] consists of an investigation of the strength of the electroweak phase transition within 
the nMSSM. In Section El we will study some of the cosmological implications of the model. In 
particular, we will investigate the constraints on the model that arise if it is to provide a realistic 
dark matter density. Section H3 consists of a discussion of the resulting phenomenology based on 
the results of the previous sections. We reserve section [7| for our conclusions. 

2 Minimal Extension of the MSSM with a Singlet Superfield 

The minimal supersymmetric extension of the SM stabilizes the gauge hierarchy, leads naturally 
to a consistent unification of the gauge couplings near the Planck scale, and provides a viable dark 
matter candidate if the neutralino, a superpartner of the neutral Higgs and gauge boson particles, 
is the LSP ^H]. Despite these successes, the model has some unattractive features as well, among 
which is the //-problem. Namely, the fi term, H2, must be included in the superpotential, with 
\fi\ of order of the weak scale, if the electroweak symmetry is to be broken. While the /i parameter is 
stable under quantum corrections as a result of supersymmetry, it is difficult (although possible 
to explain why this dimensionful quantity should be so much smaller than Mp or Mqut- 

A simple solution is to replace \i by the VEV of a gauge singlet chiral superfield, S. Other 
dimensionful couplings involving the singlet are then forbidden by demanding that the superpo- 
tential obey an additional symmetry. In the most common formulation, the Next-to-Minimal Su- 
persymmetric SM (NMSSM) 10 , one imposes a Z3 symmetry under which the fields transform as 
$j — > exp(27ri(/j/3) <3?j, where the charges % are given in tabled The most general renormalizable 
superpotential is then 

W ren = A S Hi ■ H 2 + K S 3 

+ UuQ ■ H 2 U c + ydQ ■ Hi t) c + y e L ■ Hi E c , (2) 

where S is the singlet superfield and the other fields are the same as in the MSSM. Except for 
the cubic singlet self-coupling, this is just the MSSM superpotential with a field dependent /i-term 
proportional to the singlet field. Without this additional cubic term, the superpotential is invariant 
under an anomalous U(1)pq, whose charges are listed in tabled that gives rise to an unacceptable 
axion The cubic term explicitly breaks U(1)pq down to its maximal Z3 subgroup, thereby 

removing the axion while still forbidding all dimensionful couplings. Unfortunately, this generates 
new difficulties. When the singlet acquires a VEV, necessarily near the electroweak scale, the Z3 
symmetry is broken as well producing cosmologically unacceptable domain walls |12j . The domain 
wall problem can be avoided if Z3 violating non-renormalizable operators are included. However, 
such operators generate a large singlet tadpole that destabilizes the hierarchy [T3] . 

As shown in (141 1151 ITB] . both problems can be avoided in the context of an N=l supergravity 
scenario. In the absence of the cubic singlet term, the superpotential of Eq. (j2J obeys the U(1)pq, 
and U(l)p symmetries listed in table ^ an d so is also invariant under the group U(l)^/ with 
charges R' = 3R + PQ. This symmetry alone is enough to give the superpotential of Eq. ((2j) with 
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Table 1: Charges of fields under the Abelian symmetries discussed in the text. 

no cubic term, as are the maximal Z,? and Z^ subgroups of U(l)p/. If we now demand that both 
the superpotential and the Kahler potential obey one of these discrete R-symmetries instead of 
the Z3 symmetry of the NMSSM, the superpotential will be of the desired form up to a possible 
singlet tadpole term. Using power counting arguments it may be shown that a singlet tadpole 
does arise, but only at the six (Z^) or seven (Z^) loop level. The loop suppression in both cases 
is large enough that the induced tadpole does not destabilize the hierarchy |141 I15j . Therefore, 
this mechanism very elegantly solves three problems: it prevents the appearance of dimensionful 
couplings (other than the tadpole) in the renormalizable part of the superpotential; the induced 
singlet tadpole explicitly breaks U(1)pq and its discrete subgroups, thereby avoiding unacceptable 
axions and domain walls; and the loop suppression of the tadpole leads naturally to a singlet 
VEV of the order of Msusy- Following |16| . we shall refer to this model as the nearly Minimal 
Supersymmetric Standard Model (nMSSM). 

Another interesting feature of the model is that something like R-Parity arises from the imposed 
R-symmetries. These symmetries forbid the appearance of all d = 4 B- and L-violating operators 
as well as the dominant higher dimensional S-violating operators that contribute to proton decay. 
While proton stability is ensured, there are non-renormalizable operators that make the LSP unsta- 
ble. As will be shown in Section fc.ll the LSP of the model under study is nearly always the lightest 
neutralino. In the Z^ symmetric case, this symmetry forbids all d < 6 operators that could lead 
to the decay of such an LSP, and naive dimensional analysis shows that it has a lifetime in excess 
of the age of the universe. The issue is a bit more delicate in the Z^ case since the L-violating 
d = 5 operator SSH2L is allowed by the symmetry. We find that the lifetime of the neutralino 
LSP induced by this operator is greater than the age of the Universe provided the cutoff scale 
(by which non-renormalizable operators are suppressed) exceeds A > 3 x 10 14 GeV. The details of 
our estimate are presented in Appendix lEl Therefore, the same symmetries that ensure a natural 
solution to the //-problem stabilize the LSP, and provide the means for a sufficiently large proton 
lifetime. 

The nMSSM is quite attractive for a number of phenomenological reasons as well. In the MSSM, 
large stop masses, high tan (3, and some fine-tuning are needed to evade the LEP-II Higgs mass 
bounds These bounds are much easier to avoid in the nMSSM and the NMSSM since there 
is an additional tree-level contribution to the CP-even Higgs masses. The same LEP-II bounds 
on the lightest neutral Higgs also put severe constraints on the parameter space consistent with 
Electroweak Baryogenesis (EWBG) within the MSSM 18 . On the other hand, EWBG does appear 
to be possible within the NMSSM [191 12U1 121j and other singlet extensions |22j due to additional 
terms in the tree-level potential. We find that the same holds true for the nMSSM. 

Finally, we note that the superpotential and soft-breaking terms of the nMSSM also arise as 
the low-energy effective theory of minimal supersymmetric models with dynamical electroweak 
symmetry breaking, the so-called "Fat-Higgs "models 23 . In these models, the value of the Higgs- 
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singlet coupling A is not restricted by the requirement of perturbative consistency up to the grand- 
unification scale. Instead, the precise value of A depends on the scale of dynamical symmetry 
breaking, A being larger for smaller values of this scale. In our work, we shall focus on the case 
where perturbative consistency holds up to high-energy scales, but we will also comment on how 
our results are modified if we ignore the perturbativity constraint on A. 

3 The nMSSM at Zero- Temperature 

Much of the analysis and notation in this section follows that of |15j . For simplicity, we will 
include only the Higgs, singlet, and third generation quark/squark fields in the superpotential. 
The superpotential, including the loop-generated tadpole term, is then 

2 

W nMSSM = \SH r H 2 + + y t Q-H 2 U c + . . . (3) 

where H\ = (Hi,H^\ H 2 = (h^H^j denote the two Higgs superfields, S is the singlet super- 
field, and A-B = e ab A a B b with e 12 = 1. 

The tree- level potential is then Vq = Vf + Vd + Vsoft'- 

2 

V F = \\H V H 2 + ^| 2 + \\SHl + y t t L i* R \ 2 

+ \XSH{ + y t b L t R \ 2 + \\S\ 2 h\h 2 
+ \y t f* R \ 2 HlH 2 + \y t Q-H 2 \ 2 , 

V D = ^{HlH.-HlH,) 2 + 9 ^\H\H 2 \ 2 (4) 

V soft = m\H[H x + m 2 2 H\H 2 + m 2 |S| 2 

+ {t s S + h.c.) + {a x SH r H 2 + h.c.) 

+ m 2 Q &Q + m 2 u \t R \ 2 + (a t Q-H 2 t R + h.c). 

In writing Vd we have defined g = \/ g 2 + g' 2 = gj cos 6w- 

The couplings a\,t s , \,m 2 2 ,yt,at can all be complex, but not all their phases are physical. By 
suitable redefinitions of S,Hi, and H 2 , the parameters A and m 2 2 can both be made real |24j . 
To simplify the analysis and to avoid spontaneous CP violation, we shall assume that the soft- 
breaking parameters a\,t s ,at and the Yukawa yt are also real. 1 Moreover, we may take a\ and t s 
to be positive provided we allow A and m 2 2 to have either sign. Real parameters are not sufficient 
to exclude spontaneous CP violation, however, so we must check this explicitly. We must also verify 
that the potential does not generate a VEV for either of the charged Higgs fields. 

If none of the squark fields get VEV's, the tree-level Higgs potential becomes 

V = m\H\H x + m\H\H 2 + m 2 s \S\ 2 + X 2 \H r H 2 \ 2 

lr This assumption is not completely ad hoc. Within a minimal supergravity scenario, the soft breaking parameters 
are proportional to the corresponding terms in the superpotential. 
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+ \ 2 \s\ 2 (h\h 1 + h\h 2 ) + ^{hIh 2 -h\h x ) 2 + £-\h\h 2 \ 

O Z 

+ i s (5 + /i.c.) + a x {SH v H 2 + h.c.) + m 2 12 (H 1 -H 2 + h.c). 



(5) 



We may choose an SU(2)xU(l) gauge such that (H^) = and (-Hf 1 ) E M- at the minimum of 
the potential. Taking the derivative of Vq with respect to H 2 and evaluating the result at the 
minimum, we find 



dV 



dm 



Hi=Vi 



2 i \2| |2 , 9 n |2 



2 , 5 2 , 



(6) 



where we have defined \H 2 , 



v 2 , (H£) = v+, and (5) = v s 
v\) + t vf > 0. 



It follows that (Ht) vanishes at the 



minimum provided m 2 + A 2 |v s | 2 + ^-(1^2 1 < i i > — 

If the charged Higgs VEV's vanish at the minimum, the only part of the potential that depends 
on the phases of the Higgs fields are the last three terms in Eq. ©: 



V P hase = t s (S + h.c.) + a x (SH° H Q 2 + h.c.) + mj 2 {H° H Q 2 + h.c). 



(7) 



Recalling that a\ and t s are both real and positive, the potential will have an absolute minimum 



with (S) 



v s e 



& and 



v 2 G M- provided m 2 2 < 0. While this condition is sufficient to 



avoid spontaneous CP violation, the result of indicates that it is not necessary. We will focus 
on the m 2 2 < case because it simplifies the analysis, and as we shall see below, is preferred by the 
constraints on the scalar Higgs masses. However, we have also examined the m 2 2 > case, and find 
that once we impose the experimental constraints described in the following section, the parameter 
space with m 2 2 > is very restricted, and tends to be inconsistent with electroweak baryogenesis. 
With 

m i2 ^ 0> ^ e field VEV's are all real and have fixed sign: 



(S) = v s <0 
We define the angle (3 as in the MSSM: 

v% = v cos f3, 



(Hi) = vi > 0, (tf 2 ) = v 2 > 0. 



v 2 = usin/3, 



(8) 



(9) 



with v ~ 174 GeV. We also define \x = —Xv s , since this is the quantity that corresponds to the [i 
parameter in the MSSM. Note that \x can have either sign, depending on the sign of A. 

The minimization conditions for iJP, H®, and S can be used to relate the scalar soft masses to 
the other parameters in terms of the VEV's. These give 

^2 
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-(m 12 + a x v s 
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1 dAV 



2v 2 8Hl 
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(10) 



S=v s 



where AV consists of contributions to the effective potential beyond tree-level. To one-loop order 



AV 



(4vr) s 



g b h(ml) 



\9jK 



(11) 
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where the first sum runs over all bosons, the second over all Weyl fermions, gi is the number of 
(on-shell) degrees of freedom, m^ is the field-dependent mass eigenvalue, and the function h(m 2 ) is 
given by (in the DR scheme) 



m 



him 2 ) = — 
4 

The onedoop corrections are therefore given by 
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(12) 



(13) 



Hi=Vi 



3.1 Charginos and Neutralinos 

The chargino and neutralino sectors provide important phenomenological constraints on the model. 
The fermion component of the singlet superfield, the singlino, leads to a fifth neutralino state. 
Assuming the sfermions to be heavy, with masses of order a few hundred GeV, and values of A 
that remain perturbative up to a grand unification scale of order 10 16 GeV, the LSP of the model 
is always the lightest neutralino with a mass below about 60 GeV. 
The chargino mass matrix, in the basis (W + , R^i W~ , ^T)> * s 



M x± = 





X 



X 1 




where 



X = 



M 2 V2s M w 

For the neutralinos, the mass matrix in basis ip® = (B°, W°, H®, H®, S) reads 



(14) 



(15) 



I 



\ 



Mi 


-C/3S W M Z 
Sf3S w M Z 





M 2 

C/3C W M Z 
-Sf3C w M Z 







\v 2 







\ 



/ 



(16) 



In our analysis we take M\ = ~ ^M 2 , which corresponds to what would be expected from 

universality at the GUT scale. With an eye towards electroweak baryogenesis, we allow the gaugino 
masses to have a common phase: M 2 = M e 1 ^ with M real. This phase also has a significant effect 
on the mass of the lightest neutralino. Since flipping the sign of A is equivalent to shifting the 
gaugino phase by 4> — > <p + n, we will consider only the A > case. 

To see how the light neutralino state arises, suppose M\ and M 2 are very large and real so 
that the gaugino states decouple, leaving only the lower 3x3 Higgsino block. For v\ <C v 2 ,v s , the 
smallest eigenvalue of this matrix is then 

mffi ~ 2Xv 1 v 2 v 8 /(vj + v\ + v 2 s ), (17) 

and the corresponding state is predominantly singlino. More generally, the mass eigenstates are 



Ni 



N tj 



(18) 
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where Nij is a unitary matrix such that N*Aij^N^ is diagonal with non-negative entries [I]. We 
label the states in order of increasing mass so that N\ is the lightest neutralino. 

Measurements made at LEP II impose stringent constraints on the chargino/neutralino spec- 
trum. Since the coupling of the charginos to gauge bosons is the same as in the MSSM, the mass 
of the lightest chargino must satisfy m ± > 104 GeV 1261 . The corresponding requirement for 

the neutralinos is either (m^ +m^) > 209 GeV, or a(e + e~ — > N1N2) ^ 10~ 2 fb. Finally, for 
< Mz/2, we must have BR(Z iViiVi) < 0.8 x 10" 3 

It is possible to satisfy all of these constraints in the limit of large tan (3, in which case N\ is a 
very light LSP; < 15 GeV for A < 1.0, tan (3 > 10, and M2 — > 00. This state is mostly singlino, 
and couples only weakly to the gauge bosons. However, this limit also leads to an unacceptably 
large neutralino relic density. As we will show in Section for heavy sfermions, the dominant 
annihilation channel for Ni is s-channel Z-exchange 2 . For such a light, mostly singlino JVi, the 
ZN\N\ coupling is too weak for this state to annihilate efficiently. 

We are thus led to consider values of tan/3 of order unity. The N\ state now has a sizeable 
Higgsino component and correspondingly large couplings to the gauge bosons, so there is a danger 
of producing too large a contribution to the Z-width. The branching ratio of the Z to two N\s is 
given by 

Bxz~N&-"y»r-\srfrM S < - 213/2 



4vr 2Acos 2 w T z 



1 



2m N 1 
M z 



(19) 



Combining the branching ratio constraint with that for the relic density, we find mjj ^35 GeV 
is needed if both conditions are to be met. (See Section [21) As this value depends somewhat on 
parameters in the Higgs sector, we impose the weaker constraint > 25 GeV in our analysis. 

The magnitude of A must be fairly large, A > 0.3, to raise the mass of the lightest neutralino 
above 25 GeV. (|A| > 0.5 for > 35 GeV.) On the other hand, if A is too large it encounters 
a Landau pole before the GUT scale. This is precisely what happens in the recently proposed Fat 
Higgs model |23j . in which the Landau pole corresponds to the Higgs compositeness scale. We would 
like to maintain the property of perturbative unification in the model (in the usual sense), so we 
will focus most of our attention on values of A that remain perturbative up to Mqut ~ 10 16 GeV. 
However, with the Fat Higgs model in mind, we will also consider larger values of A. 

It is straightforward to derive the limit on A at one-loop order. The relevant (one-loop) beta 
functions are ^2] 

dg s 1 3 3 

~dt = ~J4Fy2 9s, 

d Vt Vt (o 2 . 1 \2 8 2n 

-dt = W (3yt + 2 A "3 5j ' 

where t = ln(Q 2 /M^). Running these up to Q 2 = (10 16 GeV) 2 and demanding A 2 < 4-7T, we find 
the allowed region in the tan j3-X plane shown in figure ^ The lower limit on tan (3 comes about 

2 There are also contributions to the annihilation cross-section from s-channel Higgs exchange, but these processes 
alone are not strong enough to produce an acceptable neutralino relic density unless the neutralino mass is close to 
a half of one of the Higgs boson masses. 
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because small values of this quantity imply a large yt(mt), and this accelerates the running of A. 
The figure also shows the region in which the iVj state has a mass greater than 25 GeV and satisfies 
the LEP II constraints listed above. 

Figure |2 shows the corresponding allowed region in the |/u| — \M.2\ plane. The lower bounds on 
I-M2I and |/u| are due to the chargino mass constraint. Interestingly, there is also an upper bound 
on |//|, which comes from the lower bound on the lightest neutralino mass. From Eqs. (|16I17I) we 
see that for \v s \ 3> v, the predominantly singlino state becomes very light. Since [i = —Xv s and A is 
bounded above by the perturbativity constraint, this translates into an upper bound on \fi\. Both 
the phase and the magnitude of the gaugino mass M2 have a significant impact on the mass of the 
light singlino state. The largest masses are obtained for <p = tt with \M2\ ~ \v, as this maximizes 
the constructive interference between the gaugino and Higgsino components. When = the 
interference is destructive, and IM2I — * 00 maximizes the mass. 




Figure 1: Allowed regions in the tan (3— A plane. The region consistent with perturbative unification 
lies below the thick solid line, while the regions consistent with the LEP II constraints (and > 
25 GeV) lie above the thinner lines. Among these, the solid line corresponds to a gaugino phase of 
<j> = tt, while the dotted and dashed lines correspond to (j) = 0, tt/2 respectively. 

Figure El shows the range of masses of the lightest neutralino that are consistent with the 
constraints listed above. (The relevant input parameter sets are those listed in table EJ). For 
tan (3 and A below the perturbative bound we see that the N\ state has a mass below about 
60 GeV, making it the LSP in the absence of a light gravitino. For below 50 GeV, this state is 
predominantly singlino, with a sizeable Higgsino component. This is because, with the assumption 
of gaugino mass universality, the constraint on the chargino mass puts a lower bound on \M\\ that 
excludes a lighter predominantly Bino state. However, a mostly Bino LSP is possible if the light 
singlino state has a mass above about 50 GeV, although the parameter space in which this can 
occur, consistent with perturbative unification, is severely restricted. In this event, the LSP and 
NLSP must be very close in mass. If A is allowed to exceed the perturbativity bound, the situation 
is much less constrained; the parameter space in which the LSP is mostly Bino becomes much 
larger, and a Bino LSP need no longer be nearly degenerate with the NLSP. 



S 




Figure 3: Mass of the lightest neutralino. The region to the right of the solid line is consistent with 
the LEP II constraints listed above. The region surrounded by the dotted line is also consistent 
with perturbative unification. 



3.2 Higgs Spectrum 

The LEP II lower bound on the mass of the lightest neutral CP-even Higgs boson of about 114 GeV 
is difficult to evade in the MSSM. This follows from the fact that, at tree-level, the mass of this 
state is bounded by Mz, 

m 2 h < M| cos 2 2(3, (21) 

which implies that large one- loop corrections are needed to increase the mass. The dominant loop 
contribution comes from the stops. With tan/? 3> 1, stop masses of order 1 TeV, and considerable 
fine-tuning of the stop mixing parameters, the mass of the lightest Higgs can be brought up to 
m h ~ 130 GeV Ql]. 

The corresponding tree-level bound in the nMSSM is ^H] 

2A 2 

m 2 h < M\ (cos 2 2(3+— sin 2 2(3) , (22) 
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which exceeds 100 GeV for |A| ~ 0.7 and tan/3 ~ 2. The same bound applies in the NMSSM f2S] . 
This makes it possible to avoid the LEP II constraint without fine-tuning in the stop sector. 

In order to discuss some of the constraints on the parameter space from the Higgs sector, we 
list here the tree- level Higgs masses. We have also included the one-loop mass corrections from the 
top and the stops given in ^Hj in our numerical analysis. 

Since the tree-level Higgs VEV's are real, and neglecting the small CP-violating effects associ- 
ated with the one-loop chargino and neutralino contributions |29j . the Higgs fields can be expanded 
as 



Hi 
H 2 



Vl + ^(01 + Mil) 



V2 + 73(^2 + ia 2 ) 



S = v s + —7={<t>s + ia s ). (23) 
v2 

After electroweak symmetry breaking, the real and imaginary parts of the singlet mix with those 
of H® and H® to produce two neutral scalar states in addition to those of the MSSM. In all, the 
physical Higgs states consist of one charged scalar, two neutral CP-odd scalars, and three neutral 
CP-even scalars. 

For the CP-odd scalars, the combination G° = — a\ eg + a 2 sg is absorbed by the Z° while the 
orthogonal linear combination A = a\ sg + a 2 eg mixes with a s to give two physical scalars. The 
mass matrix in basis (A°,a s ) is 

Ml = ( Ml 1 u ^ 2 , ) , (24) 
p \ -a x v -l-(t s + spc B a x v 2 ) J ' v ; 

where 

M 2 a = (m? 2 + a x v s ). (25) 

S/3C/3 

Note that Eq. (JTU|) implies m? s + A 2 v 2 = —^-(t s + sg eg a\v 2 ). Therefore the singlet soft mass, m 2 , 
sets the mass scale of a predominantly singlet state. 

Among the charged Higgs bosons, the combination G + = <j>2 sg — <fy[* cp is taken up by the W + 
leaving behind a single complex charged scalar mass eigenstate, H + , of mass 

M± = Ml + Myy — X 2 v 2 . (26) 

It may be shown using the minimization conditions, Eq (|10|) . that > is equivalent to the 
condition needed to avoid a charged Higgs VEV given below Eq. ©. 

Finally, the mass matrix elements for the CP-even Higgs boson states are 

Ml = s 2 M 2 a +c 2 M 2 

Mf 2 = -sgCg(Ml + M 2 -2X 2 v 2 ) 

M 13 = v(Sf3a\ + 2cg\ 2 V s ) 

M 2 22 = 4Ml + s 2 3 M 2 

M 23 = v ( c /3 a A + ^SgX 2 V s ) 

M 3 2 3 = --(t s + s 8 cgaxv 2 )=m 2 s + X 2 v 2 (27) 

lift 
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with the remaining elements related to these by symmetry. As for the CP-odd case, the mass of a 
mostly singlet state is determined by the singlet soft mass. 

Large values of M 2 help to increase the mass of the Higgs states. This is most easily obtained 
with mf 2 < (see Eq. ([25)0 . which is a sufficient condition to guarantee the absence of spontaneous 
CP violation at tree-level. We also note that the MSSM limit of the NMSSM is not possible in 
this model. In this limit one takes \v s \ 3> v while holding \v s fixed, thereby decoupling the singlet 
states from the rest of the Higgs spectrum. As discussed in Section ETT1 such large values of v s lead 
to an unacceptably light neutralino state. On the other hand, the decoupling limit of the nMSSM 
discussed in 15], \t s \ — > oo, is still viable. Indeed, the upper bound on the lightest neutral Higgs 
mass, Eq. (|22|) is saturated in this limit if M 2 — ► oo as well. 

The precise LEP II bounds on the Higgs masses depend on the couplings of the Higgs bosons 
to the gauge bosons. These couplings tend to be weakened somewhat from mixing with the singlet. 
Let O , O p be the orthogonal mixing matrices relating the gauge and mass eigenstates: 

(tH(:;) ; (s)-^)- 

We label the mass eigenstates in order of increasing mass, so that S\ is the lightest CP-even state 
and P\ the lightest CP-odd state. In terms of these mixing matrices, the SVV-ty^e couplings are 

SZZ : 9 -M z (c p O s kl + s p O s k2 ) (Z^f S k , 

SWW : 9 -M w ( C/3 O s kl + s p O s k2 ) (W^) 2 S k . (29) 

Also relevant are the SPZ-type couplings 

SPZ: {[OUspO^-cpO^)] Z»S k d»Pt. (30) 

The couplings of the Higgs states to fermions and neutralinos are listed in Appendix O 

The LEP bound on the charged Higgs boson is given in Ref. [SO]. Assuming BR{H + — > t + u) ~ 1 
this bound reads, approximately, 

M H ± > 90 GeV. (31) 

The bounds on the CP-odd Higgs bosons depend strongly on their coupling to the Z-gauge 
boson and the CP-even scalars given in Eq. (|3U|) . If the lightest CP-odd Higgs boson, Pi, has a 
large singlet component this coupling can become very small, and the bound on this particle is 
much weaker than the LEP bound of about 90 GeV present in the MSSM [3^. This bound may 
be further weakened if the decay Pi — ► N\N\ is allowed kinematically. If so, this mode tends to 
dominate the decay width leading to a large fraction of invisible final states. 

The same is true of the lightest CP-even state, S±. The limit found in depends on the 
strength of the SVV coupling relative to the corresponding Standard Model coupling. From Eq. (|29|) 
this relative factor is {cpO^ +spOf 2 \, which can be considerably smaller than unity if the S\ state 
has a large singlet component. Again, the limit is further weakened if the S\ — > N\N\ channel is 
open, as this tends to dominate the decay width below the gauge boson threshold. In this case, the 
limit on invisible decay modes found in [HHj is the relevant one. 
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4 Elect roweak Baryogenesis 



If Electroweak Baryogenesis (EWBG) is to generate the presently observed baryon asymmetry, the 
electroweak phase transition must be strongly first order. In the most promising MSSM scenario, 
this phase transition is dominated by a light, mostly right-handed stop [HI]- Such a stop produces 
a large contribution to the cubic term in the one-loop effective potential that is responsible for 
making the phase transition first order. Even so, for Higgs masses above the LEP II bound, there 
is only a very small region of parameter space in which the EW phase transition is strong enough 
for EWBG to work QH3S]. 

The prospects for EWBG in the NMSSM are more promising. The NMSSM has an additional 
tree-level contribution to the cubic term of the effective potential. This is sufficient to make the 
electroweak phase transition strongly first order without relying on the contribution of a light stop 
|19l I2UI [2~T1 122j . Since the nMSSM has a similar cubic term in the tree- level potential, we expect 
EWBG to be possible in this model as well. 

In this section we investigate the strength to the electroweak phase transition in the nMSSM 
in order to find out whether a strongly first order transition is possible, and if so, try to map 
out the relevant region of parameter space. To simplify our analysis, we neglect the contributions 
from sfermions other than the stops since these are generally very small. We also fix the stop 
SUSY-breaking parameters to be 



This choice of parameters leads to stops that are too heavy to have a relevant impact on the strength 
of the first order phase transition. We have made this choice because it allows us to emphasize 
the effects induced by terms in the tree-level effective potential that are not present in the MSSM. 
These effects turn out to be sufficient to make the phase transition strongly first order, even in the 
absence of light stops. 

While a strongly first order electroweak phase transition is necessary for EWBG to generate 
the observed baryon asymmetry, based on previous analyses of the MSSM, it appears that this 
condition is also sufficient |36[ I37j. In the MSSM, the generation of baryon number proceeds from 
the CP-violating interactions of the charginos with the Higgs field. The dominant source of CP 
violation, leading to the baryon asymmetry, is proportional to the relative phase of the fi and the 
gaugino mass parameters, arg^fiM^). The only difference in the model under study is that the fi 
parameter is replaced by the quantity (— Xv s ), which is real, and CP violation is induced by the 
phase of the gaugino masses. Therefore, in the presence of a sufficiently strong first order phase 
transition, we expect a result for the baryon asymmetry generated from the chargino currents 
similar to the one obtained in the MSSM. 

4.1 One-Loop Effective Potential 

The finite temperature effective potential for the real Higgs scalars is 



where V n is the n-loop contribution, and the additional term, Vdaisy, wm be discussed below. Also, 
tpi, i = 1, 2, s are the classical field variables corresponding to H®, H®, and S. (jpi = V{ at the T=0 




2 Qs = m 2 % = (500 GeV) 2 
at = 100 GeV. 



V(ipi,T) = V (tpi) + Vi(ipi,T) + V daiS y(ipi,T) + ... 



(32) 
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minimum.) The tree- level part is 



Vo(ipi) = m\ipi 2 + m|</?2 2 + m 2 s ip s 2 + 2m\ 2 (pnp 2 + 2t s ip s + 2a x f s (fx (f2 



2\2 



, \2 2/ 2 i 2\ , \2 2 2 



(33) 



Note the cubic term, a\ ip s tpi <^2, which has no counterpart in the MSSM. In the DR scheme, the 
one-loop contribution reads 



V^T) = Y,9bfB(mlT)+Y,9ffF(m},T), 

b f 



(34) 



where b runs over bosons, / runs over Weyl fermions, and gi is the number of (on-shell) degrees of 
freedom. To a very good approximation, the functions fs, Jf are given by [3S] 



fB(m 2 ,T) 



fF(m 2 ,T) 



-gT* + ±m 2 T 2 - ^{m 2 ?/ 2 T - ^ H ^ 



,« 4 

64^" 



7- 



ln(2 



;m/T< 2.2 

(_^)3/2 r 4 e - m /T (1 + U T/m + _ _ ^ . m/T > 2 2 



(35) 



720 J ^ 48'" J ^ 



64^" ln (lFT 7 ' 1 
(_^)3/2 T 4 e - m /T (1 + ^ T/m + . . .) ; m /T > 1.9 



; m/T < 1.9 



647T i 



In(i 



Here, as = (47re~ 7B ) 2 , aj7 = (ire~' yE ) 2 , m 2 is the field-dependent mass at zero temperature. We 
neglect V2 and terms of higher order. 

The third contribution to the potential, Vdaisy, is a finite-temperature effect |39j . At zero 
temperature, one-loop boson self-energy diagrams which are quadratically divergent in the UV also 
develop an IR divergence as the boson mass is taken to zero. At finite-temperature, when m <C T 
the would-be IR divergence is cut off by T rather than m leading to a thermal contribution to the 
effective mass, m 2 — ► m 2 = m 2 + aT 2 . Resumming the leading thermal contributions to bosonic 
self-energies modifies the effective potential by the amount 

Vdaisy = -Tg- 9b ~ m b) 3/2 > ( 36 ) 



where m 2 is the thermal mass. The sum includes gauge bosons, although only the longitudinal 
modes of these develop a thermal contribution to their mass at leading order. The field-dependent 
mass matrices relevant to our analysis (including thermal corrections) are listed in Appendix 1X1 



4.2 Tree-Level Analysis 

To better understand the effect of the new cubic term, we have examined a simplified form of the 
potential which allows us to obtain analytic expressions for the critical temperature, T c , and the 
field VEV's. If the cubic term plays a dominant role in making the electroweak phase transition 
first order, we expect this analysis to give a good qualitative description of the transition. 

Our first simplifying assumption is that the ratio of the field values at the broken phase minimum 
remains constant up to T c . That is, we fix tan /3, and consider variations in ip = \J (pi 2 + ip2 2 - To 
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make the one-loop part of the potential more manageable, we keep only the leading ip 2 T 2 terms 
in the low-temperature expansion, and include only the contributions of the gauge bosons and the 
top. Since the stops are assumed to be heavy, the leading temperature-dependent contribution 
comes from the top provided A lies below the perturbative bound, A < 0.8. For larger values of A 
the contributions from the Higgs, charginos, and neutralinos become important. We shall restrict 
ourselves to the perturbative regime in the present analysis. 

In terms of if, <p s , and /3, the effective potential becomes V = Vo + V\. The tree- level part is 

V = MV + m 2 s <Ps 2 + 2t s (p s + 2a<p 2 (p s + \ 2 <p 2 <p s 2 + A V, (37) 
where we have defined 

M 2 = m\ cos 2 (3 + ml sin 2 (3 + m\ 2 sin 2/3, 
a = a\ sin (5 cos (3, 
X 2 a 2 

A 2 = — sin 2 2/3 + — cos 2 2/3. (38) 
4 8 

Within our approximation, the one-loop part is 



Vi = 1[9 2 + y + 2y 2 sin 2 /3j^ 2 T 2 . (39) 

To find the minimization conditions, we shall make use of the simple (quadratic) cp s dependence 
and consider only the field-space trajectory = along which the minimum of the potential is 
found. This condition allows us to eliminate ip s in terms of ip giving 



-,2 



^m 2 + \ 2 ip 2 

Inserting this back into the effective potential, we find 

V(<p, T) = m W " + XV ' ( 41 ) 

where m 2 (T) = M 2 + ± (g 2 + ^ + 2y 2 sin 2 ^ . 

The critical temperature, T c , and the VEV at T c , ip c , are defined by the two conditions 

V(^,T C ) = V(<p = 0,T c ), (42) 
dV 



dip 



0. 



*P=<Pc 



Solving for ip c and T c we find 



(p 2 . = —r[-m^ + —\m s d 



If 9 1, ~ X 2 t s 



A 2 V s A m s 

(43) 



T 2 = 8(F( 9 , 2 )-F(, 2 ))/( 5 2 + | + 2 y 2 sin 2 /3), 
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where 



2d 



t s + cup* 
m 2 + \ 2 ip 2 



t s + d(p 2 
m 2 + A 2 (/9 2 



2,„2 



2AV 



(44) 



Both </? 2 and T 2 must of course be positive if a solution is to exist. For <p 2 > 0, we need 



2 1 

m„ < — 

A 



A 2 i.< 



m s a 



(45) 



The positivity of T 2 requires that F(ip 2 ) > F{v 2 ). Since increasing the temperature tends to 
decrease the field VEV, this condition will be satisfied if F(ip 2 ) is a decreasing function which is 
the case provided 

X2 ^-) 2 <A 2 (m 2 + AV) 3 . (46) 



[m s a 



m. 



It is sufficient to demand that this hold for ip = ip c , which gives 

m 3 (m 2 d - X 2 t s ) 2 < \\m 2 s a - A 2 t s | 3 . 
A 



(47) 



This is equivalent to the inequality in Eq. (|45|) . Thus, Eq. I)45|) is the necessary condition for a first 
order phase transition. To satisfy this equation, m 2 must not be too large. This can be also seen 
from Eq. (|41j) . which has only positive quadratic and quartic terms in the limit A 2 ^ 2 <C m 2 . 



4.3 Numerical Analysis 

The results of the previous section have been examined more carefully by means of a numerical 
investigation of the one-loop effective potential. In this analysis, we consider only the dominant con- 
tributions which are those of the top, the stops, the gauge bosons, the Higgs bosons, the charginos, 
and the neutralinos. The corresponding field-dependent mass matrices, both at zero and finite 
temperature, are listed in Appendix EI For the purpose of calculating thermal masses, we assume 
that the remaining sfermions and the gluino are heavy enough to be neglected. We find that a 
strongly first order electroweak phase transition is possible within the nMSSM. 

The procedure we use goes as follows. To begin, we specify the values of 
((3, v s , a\, t s , M a , A, I-M2I, (/>), where M2 = \M2\e % ^ is the complex Wino soft mass. These are 
chosen randomly from the initial ranges listed in table [21 As above, the Bino mass is taken to 
be Mi = M2/2, and we fix the soft stop parameters to be TOq 3 = m 2 ^ = (500 GeV) 2 , at = 
100 GeV. The subtraction scale is set at Q 2 = (150 GeV) 2 . For each parameter set we calculate 
the mass spectrum at zero temperature and impose the experimental constraints described in 
Section |21 At this point we do not impose any dark matter constraint other than the necessary 
condition > 25 GeV. Since, for some parameter sets, the one- loop correction can destabilize 
the potential in the ip s direction, we also check that the minimum at (cpi, (p%, (p s ) = (vi,V2,v s ) is 
a global minimum at T = 0. Finally, we calculate (p c and T c using the full potential, where T c is 
taken to be the temperature at which the symmetric and broken phase minima are equal in depth, 
and ip c = \J (fi 2 + (f2 2 is the broken phase VEV at this temperature. 

In this way we have found several parameter sets which satisfy all the constraints listed above, 
and give 9? C /T C > 0.9, which we take as our criterion for a strongly first order transition |34| I40j . 3 

3 This corresponds to ^ > 1.3 for ip normalized to 246 GeV. 
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tan (5 


A 


V s 

(GeV) 


a\ 
(GeV) 


.1/3 

(GeV) 


M a 
(GeV) 


|M 2 | 
(GeV) 




Range 


1-5 


0.3-2.0 


-750-0 


0-1000 


0-1000 


200-1000 


0-1000 


0-7T 



Table 2: Ranges of input parameters. 



Let us define 



D 



X ml 



X 2 U 



a m< 



(48) 



where m 2 = — a a v\ V2/v s — t s /v s — X 2 v 2 (see Eq. (|10j)). and a and A are defined in Eq. 1)38(1 . D is 
crucial in determining whether or not a first order phase transition occurs. The simplified analysis 
of Section 14.21 Eq. (|45|) in particular, suggests that \D\ > 1 is a necessary condition for a first 
order transition. Figure |1] shows D plotted against m s for both A below the perturbative bound, 
and for general values in the range 0.7 < A < 2.0. The region surrounded by the solid lines in 
this figure corresponds to parameter sets consistent with the experimental constraints. This figure 
shows that among the parameter sets for which a strongly first order phase transition occurs, most 
satisfy \D\ > 1. On the other hand, we also find parameter sets with \D\ > 1 that do not exhibit 
a strongly first order phase transition, so this condition is not a sufficient one. The low m s region 
in these plots is excluded since the potential tends to become unstable in the singlet direction for 
small values of this quantity. This leads to the additional requirement of m 2 > (50 GeV) 2 . 



D 





Figure 4: Values of D for parameter sets leading to a strongly first order phase transition for: (a) 
A below the perturbative bound; (b) general values of A in the range 0.7 < A < 2.0. The region 
consistent with the experimental constraints lies within the area enclosed by the solid lines. 

The critical temperature for the phase transition generally falls in the range T c = 100—150 GeV. 
Table |21 shows the parameter values and transition temperatures for three of the successful param- 
eter sets. The particle spectra corresponding to these are listed in Appendix O Parameter sets A 
and B, with \D\ ~ 1, both satisfy the perturbative bound while C, for which D ~ 6.7, exceeds it. 

Sets A and B are typical of the (perturbative) parameter sets that give a strong phase tran- 
sition. As we found in Section l3.1| the constraints in the chargino/neutralino sector, along with 
perturbative consistency, force A ~ 0.5— 0.8, tan/3 ~ 1.5 — 5, and \v s \ ~ 150 — 500 GeV. For a 
given v s , the values of a\ and t s must then be adjusted so that m 2 > (50 GeV) 2 (Eq.(^J)) arid 
\D\ > 1 (Ea. H48|0 if the potential is to be stable and the transition is to be strongly first order. 
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These quantities are further constrained by the Higgs sector. We find that these requirements may 
be satisfied for a\ ~ 300—600 GeV and t s ~ (50—200 GeV) 3 . An upper bound on the value of m s 
is obtained, that is about 200 GeV. This bound arises from the phenomenological constraints and, 
most importantly, due to the dependence of the parameter D on m s , Eq. (|48|). from the condition 
\D\ > 1. This bound on m s has important consequences for Higgs physics, as we will describe in 
section El 

The value of M a , instead, does not appear to have much effect on the phase transition, but 
tends to be fairly large, M a > 400 GeV, due to the Higgs mass constraints. While large values of 
M a help to increase the masses of the lightest Higgs states, they also tend to make EWBG less 
efficient [351 1361 13*7| . Even so, EWBG is still able to account for the baryon asymmetry provided 
tan/3 < 2, as we tend to find here [35] . 



Set 


tan (3 


A 


v s 
(GeV) 


a\ 
(GeV) 


1/3 

(GeV) 


M a 
(GeV) 


\M 2 \ 
(GeV) 





(GeV) 


T 
(GeV) 


Qh 2 


A 


1.70 


0.619 


-384 


373 


157 


923 


245 


0.14 


120 


125 


0.102 


B 


1.99 


0.676 


-220 


305 


143 


914 


418 


2.57 


145 


95 


0.010 


C 


1.10 


0.920 


-276 


386 


140 


514 


462 


2.38 


145 


130 


0.104 



Table 3: Sample parameter sets exhibiting a strongly first order electroweak phase transition. 



5 Neutralino Dark Matter 

As discussed above, the LSP in this model (for A below the perturbative bound) is always the 
lightest neutralino, N%, with a mass below about 60 GeV and a sizeable singlino component. This 
can be dangerous since a light, stable particle with very weak gauge couplings may produce a relic 
density much larger than is consistent with the observed cosmology. On the other hand, if the N\ 
is able to annihilate sufficiently well, this state makes a good dark matter candidate |41j . 

For values of tan/3 and A consistent with the perturbative limit, the LSP tends to be mostly 
singlino, but has a sizeable higgsino component. Since TUft < 60 GeV, s-channel Z° exchange is 
the dominant annihilation mode. There are also contributions from s-channel CP-even and CP-odd 
Higgs boson exchanges generated by the the A S H\ ■ H 2 term in the superpotential, although these 
tend to be very small except near the corresponding mass poles. Since we have assumed that all 
sfermions are heavy, we consider only these two channels in our analysis. The relevant couplings 
are listed in Appendix El (The matrix element for s-channel Z-exchange is given in |42| .) Let us 
stress that coannihilation between the lightest neutralino and the other chargino and neutralino 
states is not important, since for the parameter sets consistent with EWBG and the perturbative 
bound, the NLSP is always at least 15 % (and almost always more than 25%) heavier than the 
LSP, implying that the coannihilation contribution is Boltzmann-suppressed. 

The LSP mass range of interest lies near the Z-pole making the annihilation cross-section 
a rapidly varying function of the mass. Since this can cause problems for the non-relativistic 
expansion commonly used to calculate the thermal average of the cross-section [13] > we have instead 
followed the methods described in [11] to do the thermal average. This gives 

roc _ . 

(av}= ds \/s - 4M 2 W K x (y/s/T) /WM 4 T K 2 (M/T) , (49) 
Jam 2 1 
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where M = is the neutralino mass, T is the temperature, s is the usual Mandelstam parameter, 
K\ and K2 are modified Bessel functions, and the quantity W is defined to be 



W 



n 



d 3 Pf 



(2n) 3 2E f 



f 



(50) 



where \Ai\ 2 is the squared matrix element averaged over initial states, and summed over final states. 

To find the relic density we have solved the corresponding Boltzmann equation using the ap- 
proximation method described in |45| . The ratio of mass to freeze out temperature, xt = M/Tf, is 
given by the solution to 



In 



0.038 {g/gl / s 2 )MM Pl (ay) (x f ) 
1/2 



(51) 



where g = 2 is the number of degrees of freedom of the neutralino, g* s is the total number of 
relativistic degrees of freedom at temperature Tf (we use a simple step approximation for both g* 
and g* s ), and Mpi is the Planck mass. The relic density is then given by 



O/i 2 



;i.07 x 10 9 GeV -1 ) 
Mpi 



x f 



(52) 
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Figure 5: Neutralino relic density as a function of mass for two values of the mixing parameter, 
7Vi3 1 2 — |iVi4| 2 | = 0.1 (dotted), 0.5 (dashed) and typical values of the Higgs mixing parameters. 
The region to the right of the thick solid line is consistent with the observed Z width. The scattered 
points correspond to parameter sets that give a strong first order phase transition, and are consistent 
with perturbative unification. 

Figures El and IH1 show the relic densities obtained for parameter sets that satisfy the abovemen- 
tioned experimental constraints, and are consistent with EWBG. A relic density consistent with 
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Figure 6: Neutralino relic density as a function of mass for two values of the mixing parameter, 
||-^13| 2 — I-^14-| 2 | = 0.1 (dotted), 0.5 (dashed) and typical values of the Higgs mixing parameters. 
The region to the right of the thick solid line is consistent with the observed Z width. The 
scattered points correspond to parameter sets which give a strong first order phase transition with 
0.7 < A < 1.2. 

the observed dark matter is obtained for neutralino masses in the range rrifj ~ 30 — 40 GeV. For 
neutralino masses greater than this, it appears to be difficult to generate a sufficiently large relic 
density to account for the dark matter. If we allow A to exceed the perturbativity bound, a realistic 
dark matter candidate may be obtained for neutralino masses above about 65 GeV. 

6 Phenomenological Discussion 

The region of parameter space consistent with EWBG, neutralino dark matter, and the experi- 
mental bounds is quite constrained, and leads to an interesting phenomenology. We shall focus on 
values of tan (5 and A that satisfy the perturbativity bound. 

The dark matter condition implies that the LSP of the model is the lightest neutralino with 
a mass in the range rrifj — 30 — 40 GeV, and is mostly singlino. For smaller values of I-M2I, the 
next-to-lightest neutralino is predominantly bino. Otherwise it is a mostly Higgsino state. In both 
cases, there are always two mostly Higgsino states with masses of order \\v s \ < 350 GeV. The 
bound on the Higgsino masses comes from the bound on found in Section f3.ll This bound also 
implies that the lightest chargino state has a mass below this value. 

In the Higgs sector, since M a tends to be fairly large, one CP-even Higgs state, one CP-odd 
Higgs state, and the charged Higgs end up with large masses of order M a . The remaining CP-odd 
state is relatively light, and is nearly a pure singlet. For M a — > 00, the tree-level mass of this state 
goes to (see (Eq. E3J)) 

m P -> ^-(t s + a x v 2 s l3 c j3 )/v s = ^ m 2 + X 2 v 2 , (53) 
which is less than 250 GeV for the values of t s , v s and a\ consistent with EWBG, as may be seen 
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from Figure |1J 

The remaining two CP-even states also tend to be fairly light, 
effective tree-level mass matrix for these states becomes 



In the M, 



oo limit, the 



Ml 



( 




| cos 2 2(3 + A 2 v 2 sin 2 2(3 v(a x sin 2(3 + 2 X 2 v s ) 
v(a x sin 2/3 + 2 X 2 v s ) -(t s + a x v 2 sp cp)/v. 



) 



(54) 



Sl,2 ~ 



s 



Among the parameter ranges consistent with EWBG, the off-diagonal element of this matrix can 
be of the same order as the diagonal elements which leads to a strong mixing between the gauge 
eigenstates. However, too much mixing can produce an unacceptably light mass eigenstate, in 
conflict with the LEP bounds [321 I33| . To avoid this, the two terms in the off-diagonal matrix 
element must cancel to some extent. If a\ = — 2X 2 v s / sin 2/?, the cancellation is exact, and the 
mixing goes to zero. The mass eigenstates then consist of the SM-Higgs-like linear combination 
(cos (3 + sin/3 H®) with a tree- level mass below 115 GeV, and a singlet state that is degenerate 
with the lightest CP-odd state. In general these states do mix, but since the mixing is limited 
by the Higgs mass constraint, one state remains predominantly SM-like while the other is mostly 
singlet. For finite M a , the corrections to this picture are of order a x v 2 /M 2 . 

The mass of the SM-like CP-even state is increased significantly by the one-loop corrections 
from the top and the stops. For tuq = mfj = 500 GeV, at = 100 GeV, as used in our analysis, we 
find that the mass of the lightest Higgs boson can be as large as 130 GeV, and still be consistent 
with EWBG. 4 This is somewhat larger than the corresponding MSSM limit of about 120 GeV |Hj, 
and is the result of the additional terms in the tree-level potential. These lead to a larger Higgs 
quartic coupling than in the MSSM, which sets the scale of the lightest SM-like CP-even state, 
while still allowing for a strongly first order electroweak phase transition. 

Figures and |H1 show the mass and composition of the three light Higgs states (including 
the one-loop contributions from the top and the stops) for the representative parameter values 
M a = 900 GeV, t s = (150 GeV) 3 , v s = -300 GeV, a x = 350 GeV, and A = 0.7. These values 
are typical of those consistent with the constraints and EWBG. The maximum of the mass of the 
lightest CP-even state occurs when the off-diagonal term in Eq. (|54|) vanishes; sin 2(3 = — 2A 2 v s /a\. 
The splitting between the Pi and S2 states at this point is due to the finite value of M a . As tan (3 
varies away from the maximum point, the mixing between the CP-even Higgs boson states increases 
and, for the specific parameters chosen in Fig. values of tan (3 > 2.5 are excluded by the current 
Higgs mass bounds obtained by the LEP experiments [3*2*1 13*3] . 

The discovery of these light states is more challenging than the MSSM case for two reasons. 
First, all three can have sizeable singlet components which reduce their couplings to the gauge 
bosons and quarks, and therefore their production cross-sections. Second, these states can decay 
invisibly into pairs of the neutralino LSP. For CP-even Higgs masses below the weak gauge boson 
production threshold this mode dominates the decay width: BR(S —* N\Ni) ~ 0.60 — 0.95 for the 
SM-like state; BR(S — > N\Ni) ~ 1 for the mostly singlet state. However, for masses above (or 
near) the threshold, weak gauge boson final states become dominant unless the Higgs boson is a 
nearly pure singlet state. 

Previous studies indicate that the most promising discovery mode for an invisibly decaying CP- 
even Higgs at the LHC is vector boson fusion (VBF) |4f)| l47 | H%] . In [IB] the authors find that, from 
the invisible modes alone, Higgs masses up to 150 GeV can be excluded at 95% CL with 10 fb _1 

4 We do not expect two-loop corrections to significantly alter this result for the set of stop parameters used here. 
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Figure 7: Mass and of the light Higgs bosons for typical parameter values consistent with EWBG. 
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Figure 8: Composition of the two light CP-even Higgs bosons for typical parameter values consistent 
with EWBG. The solid lines correspond to the components of Si, Of i: while the dotted lines are 
those of S*2, Ofj, for i = 1,2, S. The composition of the Pi state is not shown because it is almost 
pure singlet. 



of integrated luminosity provided rj > 0.35, where 

r, = BR{h^im) °%BF) (g5) 

(7(1/ Bb )SM 

If we treat these limits as being due only to statistical uncertainty and rescale them by luminosity, 
we obtain 

1.2 fb' 1 8.0 fb' 1 

£ 95% - 2 , £5* - ^ , (56) 

where £95% and £50- are the luminosities needed for a 95% CL exclusion and a 5a discovery 
respectively. For the SM-like CP-even state, we find 77 ~ 0.5 — 0.9 among the parameter sets 
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consistent with EWBG. This implies that, from the invisible channels alone, less than about 5 fb~ 
of integrated luminosity is needed to exclude this state at 95% CL, while 10 — 30 fb _1 is sufficient 
for a 5a discovery. Similarly, we find r\ ~ — 0.35 for the mostly singlet CP-even state if the mass 
lies below the gauge boson threshold. Thus, at least 10 fb _1 is needed for a 95% CL exclusion and 
65 ftT 1 for a 5a discovery using the invisible modes. On the other hand, if this state has a mass 
above the gauge boson threshold, the Higgs component is usually large enough that gauge boson 
final states completely dominate the branching ratio. Lastly, the light CP-odd state is nearly pure 
singlet and tends to decay invisibly, making it extremely difficult to observe. 

7 Conclusions 

The origin of the matter-antimatter asymmetry and the source of the dark matter are two of 
the most important questions at the interface of particle physics and cosmology. In this article, 
we have investigated these questions within a next-to-minimal super symmetric extension of the 
Standard Model. The nMSSM model, which elegantly solves the fi problem by adding a gauge 
singlet superfield, appears to be consistent with all current experimental constraints, does not 
suffer from the usual domain wall problem of the NMSSM, and leads to the stability of the proton 
and a neutralino LSP over cosmological times. 

We have shown that a strongly first order electroweak phase transition, necessary to preserve 
the baryon asymmetry produced by Electroweak Baryogenesis, may be naturally obtained within 
this model. The strength of the phase transition is largely determined by terms in the tree-level 
nMSSM scalar potential. This differs from the MSSM, in which one-loop corrections from a light 
stop are needed to make the transition first order. 

We have also shown that, if perturbative consistency is required to hold up to the GUT scale, 
the LSP of the model (in the absence of a light gravitino) is always the lightest neutralino, and 
has a mass below about 60 GeV. This state provides a viable dark matter candidate for neutralino 
masses in the range 30-40 GeV. Furthermore, we find that it is possible to achieve simultaneously 
both a realistic neutralino relic density and a strongly first order phase transition. 

In the region of parameters consistent with both requirements, there are always at least two light 
CP-even and one light CP-odd Higgs bosons. These can decay invisibly into neutralinos providing 
an interesting modification to the standard Higgs physics processes. 
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A Field-Dependent Masses 

In this section we collect the field and temperature dependent mass matrices for those particles 
relevant to the analysis in Section 0J The leading thermal mass corrections were calculated follow- 
ing U9] for vanishing background field values; (fi = cp% = ip s = 0. Ignoring the singlet background 
is reasonable here since, in the parameter space of interest, the singlet VEV is closely related to the 
other Higgs VEV's. To leading order, only bosons receive thermal mass corrections. These come 
from quadratically divergent (at T=0) loops containing particles which are light compared to the 
temperature; m < 2ttT for bosons and m < ttT for fermions. For the purpose of calculating ther- 
mal masses we have taken the Higgs bosons, Higgsinos, electroweak gauginos, and the SM particles 
to be light, while treating the rest of the particles in the spectrum as heavy. We do not expect that 
a different choice of spectrum would change our phase transition results since the first order nature 
of the transition is determined by the tree-level potential rather than the cubic one-loop term in 

3B- 

A.l Gauge Bosons 

At leading order, only the longitudinal components of the gauge bosons receive thermal corrections. 
The masses are 




(57) 



where IT = for the transverse modes, and 




6 



(58) 



for the longitudinal modes. 



A. 2 Tops and Stops 




(59) 



where 






(60) 
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A. 3 Higgs Bosons 

The CP-even mass matrix elements are 

a 2 

M 2 Sll = mj + X 2 ^ 2 + i Ps 2 ) + ^(3i Pl 2 -ip 2 2 ) + U Hl 
M 2 Sl2 = mj 2 + 2i Pl ^ 2 (X 2 - —) +a x ^ s 

-M-Sm = °A^2 + 2\ 2 (p 1 (p s 

a 2 

M% 22 = m 2 + \ 2 (cp 1 2 + cp s 2 ) + ^(3i P2 2 -i Pl 2 ) + Il H2 

M 2 S2:i = a\(p! + 2\ 2 (p 2 ip s 

M 2 s 33 = m 2 + X 2 ( Vl 2 + ^ 2 2 )+U s 

where the leading thermal corrections are 

= \ 9 ,2 t 2 + \g 2 T 2 + 1a 2 t 2 

% 2 = \g ,2 T 2 + \ 9 2 T 2 + ±\ 2 T 2 + \y 2 

lis = \>?T 2 . 
6 

For the CP-odd states we have 





Q 2 


-^ 2 ) + U Hl 




= ~ m 12 ~ a Ws 




M 2 Pl3 


= -a\tp 2 




M 2 P22 


Q 2 

= ml + X 2 ^ 2 +^ 2 ) + ^( V 9 2 2 


- ^i 2 ) + U H2 




= -a\(fi 




M 2 P33 


= m 2 + X 2 ^ 1 2 + i P2 2 )+Us 





where the thermal corrections are as above. 
The charged Higgs mass matrix is 

-2 2 

M%± = ml + xW - 9 -{^ 2 -^ 2 ) + 9 -^ 2 + n Hl 
g 2 

M 2 H ± = -(A 2 - —)(fi(p2 - m\ 2 - a x ip s 

a 2 a 2 
M H± = rn 2 + X 2 ^ s 2 + ^ 2 2 - Vl 2 ) + ^ l 2 + Ii H2 . 

A. 4 Charginos and Neutralinos 

The chargino mass matrix reads 

M * = {x X o) 



24 



where 

x = ( M 2 g(f2 \ 
V 9<Pi -^<Ps J ' 

For the neutralinos we have 



/ 


Mi 








■ \ 







M 2 










~fa 


fa 











fa 


-fa 


Xlfs 







V 








\lf2 




o / 



We have taken Mi = M 2 /2 and have allowed M 2 to be complex in our analysis. 



B Higgs and Neutralino Couplings 

We list here the Higgs and Neutralino couplings relevant to our analysis. All fermions are written 
in terms of four-component spinors to facilitate the derivation of the Feynman rules. As above, 
Eq. (|18jl. we define the unitary matrix Ny by 

Nt = Nij $ (68) 

where (ipi) = (B°,W°, H®, H®, S). Similarly, as in Eq. 1)28(1 we define the orthogonal matrices O s 
and O p by 




B.l Neutralino Couplings 

ZNN: [J 

&znn = ^-^-Z^N, -fiOZPL - OfPajNj (70) 

l cos UW 

where 

°ij = \ (^iV; 4 - N i3 N* 3 ) . (71) 

For the Higgs-Neutralino couplings we consider only the contribution from the \SH\-H2 term 
in the superpotential. We therefore neglect the contribution due to mixing with the gauginos. 
SNN: _ 

XSNN = -^ S ^i ( A ij P L + ^i' ( 72 ) 

where 

Aij = ®kl Qij + ®k2 Qij + Cfc3 Qij) 

Qt 3 = l{N: a N* b + N; b N; a ). 

A|- is related to A\- by the replacement N*- — > N^. 
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PNN 

= -i _ 
V2 



where 



&pnn = -i—S k Ni B%P L - B%P R Nj, (73) 



Pff- 



B% = spO^Qff + c /J 0&Qg i + OgQ%, 

As with Af H , Bf 4 is related to Bf- by the replacement N*- — ► Na . 

Finally we note that in converting these couplings into Feynman rules, one must insert an 
additional factor of two since the neutralinos are written as Major ana spinors pQ. 

B.2 Higgs Couplings 

The relevant couplings of the Higgs scalars to the gauge bosons are given in Section EI3 above. We 

list here the couplings of the Higgs to the third generation quarks. 

Sff: 

Se pff = ^VbSp 0& P k h 5 b +^Vtc P O s kl P k Ht. (75) 
The couplings of the Higgs states to the other matter fermions follow the same pattern. 

B.3 LSP Lifetime 

Having listed the Higgs-Neutralino couplings, we may now present an estimate for the lifetime 
of the LSP in the symmetric scenario. This symmetry allows the d = 5 operator SSH2L in 
the superpotential, which can lead to the decay of the neutralino LSP. This operator generates a 
coupling which allows the neutralino to decay into a neutrino and a pair of off-shell neutral Higgs 
scalars, or an electron, a neutral Higgs scalar, and a charged Higgs scalar. We shall focus on the 
first mode with intermediate neutral CP-even states and a mostly singlino LSP. This gives a more 
stringent constraint than the charged mode, and the analysis with intermediate neutral CP-odd 
states is analogous. We will also assume that each of the neutral Higgs bosons subsequently decays 
into bb. With these assumptions, we find 



7T 



Vb 



4 at* f\S snS snS snS 

l '" Vl - """"" ' 2<w N 1 ) 



m<- T \ 11 



2 



where A is the cutoff scale, ran is the Higgs mass in the intermediate propagators, and we have 
set all final state momenta to Setting the mixing factor to unity, taking m^/mH ~ 1 and 

tan/3 = 2, and demanding that V < Hq, we find 

A 2>^) 3 1 23 GeV 2 (77) 



GeV 

which translates into A > 3 x 10 14 GeV for = 100 GeV. 
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C Sample Mass Spectra 

We list here the mass spectra for the sample parameter sets A, B, and C listed in table 01 
Higgs Scalar Masses: 



Set 


Si 


s 2 


S3 


Pi 


P 2 


H± 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


A 


115 


158 


925 


135 


927 


922 


B 


116 


182 


914 


164 


917 


911 


C 


121 


219 


504 


115 


534 


498 



Neutralino and Chargino Masses: 



Set 


Ni 


N 2 


^3 


N 4 


N 5 


xt 


xt 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


A 


33.3 


107 


181 


278 


324 


165 


320 


B 


52.4 


168 


203 


221 


432 


151 


432 


C 


77.1 


228 


268 


331 


474 


257 


474 



For completeness, we also display the composition of the lightest neutralino, N±, the lightest 
CP-even Higgs, Si, and the lightest CP-odd Higgs, Pi. 



Set 




\N l2 \ 


\N 13 1 


\N U \ 


1^15 1 


(OS 


u 12 


^13 


^11 


O p 

u 12 


A 


0.13 


0.10 


0.11 


0.37 


0.91 


0.46 


0.74 


-0.50 


0.08 


0.99 


B 


0.07 


0.07 


0.16 


0.52 


0.84 


0.42 


0.80 


-0.44 


0.07 


0.99 


C 


0.01 


0.01 


0.28 


0.33 


0.90 


0.70 


0.53 


0.48 


0.26 


0.97 
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